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We introduce a simple yet powerful analytic method which obtains the structure of cosmic microwave back- 
ground anisotropics to better than 5-10% in temperature fluctuations on all scales. It is applicable to any 
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^^ model in which the potential fluctuations at recombination are both linear and known. Moreover, it recov- 

(^— ^ ers and explains the presence of the "Doppler peaks" at degree scales as driven acoustic oscillations of the 

^^ photon-baryon fluid. We treat in detail such subtleties as the time dependence of the gravitational driving 

"^ force, anisotropic stress from the neutrino quadrupole, and damping during the recombination process, again 

O I all from an analytic standpoint. We apply this formalism to the standard cold dark matter model to gain 

I 

Q physical insight into the anisotropics, including the dependence of the peak locations and heights on cosmo- 

5— I . 

-^-^ logical parameters such as fi;, and h, as well as model parameters such as the ionization history. Damping 

C^ due to the finite thickness of the last scattering surface and photon diffusion are further more shown to be 

identical. In addition to being a powerful probe into the nature of anisotropics, this treatment can be used 

in place of the standard Boltzmann code where 5-10% accuracy in temperature fluctuations is satisfactory 

and/or speed is essential. Equally importantly, it can be used as a portable standard by which numerical 

codes can be tested and compared. 
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It is the nature of things that they are ties to each other. 

-Chuang-tzu 

1. Introduction 

After their discovery by the COBE satellite (Smoot et al. 1992), cosmic microwave background (CMB) 
anisotropics have become one of the most powerful observational probes of cosmology. Indeed recent detec- 
tions of CMB anisotropics on degree scales (see e.g. White, Scott, & Silk 1994 and references therein) provide 
us with important information about large scale structure formation in the universe. It is consequently of 
great interest to understand their origin. To predict CMB anisotropics in any given model, we have to solve 
the coupled equations for the evolution of all species present. Namely, these are the Euler and continuity 
equations for the fluid components (the baryons and cold dark matter) and the Boltzmann equations for the 
massless components (the photons and neutrinos). These coupled equations have been numerically solved 
by many authors {e.g. Peebles k Yu 1970, Wilson k Silk 1981, Bond k Efstathiou 1984, Vittorio k Silk 
1984). Although it is sufficient for direct comparison of a specific model with observations, this "black box" 
approach makes it difficult to extract the physical content of the predictions. This problem is compounded 
by the fact anisotropy formation is a rather complicated process. It has long been known that several phys- 
ically distinct effects contribute to their generation, e.g. the gravitational redshift (Sachs k Wolfe 1967), 
the adiabatic growth of perturbations (Peebles k Yu 1970), baryon velocity induced fluctuations (Zel'dovich 
k Sunyaev 1970), and photon diffusion (Silk 1968). An analytic treatment is therefore desirable to gain 
physical insight into CMB anisotropics. 

In this paper, we present a fully analytic treatment for the evolution of CMB temperature perturbations 
and systematically investigate each contribution to the final observable anisotropy. It can be applied to any 
model with the standard thermal history, regardless of dark matter content or initial conditions. Several 
analytic calculations of CMB anisotropics have been performed in the past under less general, and often 
unrealistic, assumptions. Doroshkevich, Zel'dovich, k Sunyaev (1978) and Doroshkevich (1988) presented 
an analytic expression for temperature anisotropics on the last scattering surface. Based on this work, 
Naselsky and Novikov (1993), J0rgensen et al. (1994), and Atrio-Barandela k Doroshkevich (1994) have 
recently calculated the anisotropics in the cold dark matter (CDM) scenario. However these works did not 
realistically account for the evolution of the gravitational potential inside the Jeans scale during radiation 
domination or on any scale through matter-radiation equality. Indeed CMB anisotropics are quite sensitive 
to this evolution as we shall show. Thus agreement with the numerical solutions could not be established 
in those works. Furthermore, they used an overly simplistic account of fluctuation evolution during the 
recombination process when the damping scale of CMB anisotropics is fixed. They also neglected the 
neutrino contributions to the anisotropic stress which is important during radiation domination. 

Moreover the synchronous gauge condition employed by most previous treatments makes the separa- 
tion and physical interpretation of these effects difficult. Here we use gauge invariant perturbation theory 
(Bardeen 1980; Kodama k Sasaki 1984) where each physical process is readily distinguished. Recently a 
gauge invariant treatment of CMB anisotropics has been performed by Seljak (1994) based on similar approx- 
imations to our own. However, it employs a numerical solution to the tightly coupled evolution equations. 
Our analytic treatment allows one to separate and individually examine each contribution to the anisotropy 
easily, yet still maintains sufficient generality to realistically describe a model such as CDM. 

In §2, we introduce the general technique and present the central results of our analytic approxima- 
tion, valid for any gravitational instability model. Starting from the exact multifluid differential equations 
(§2.1), we derive the tight coupling approximation (§2.2) and obtain its analytic solution in terms of simple 
quadratures (§2.3). The tight coupling approximation is modified to include photon diffusion through re- 
combination in §2.4. All of these approximations are treated in greater detail in Appendix B. Furthermore, 



useful analytic formulae which describe the recombination process are presented in Appendix C. In §2.5, we 
discuss the free streaming solution and show how the final observable anisotropics are obtained. 

We specialize these techniques to the CDM model in §3. The so-called "Doppler peaks" in the anisotropy 
spectrum are shown to arise from driven acoustic oscillations in the tight coupling regime. The analytic form 
of the gravitational driving force is derived in Appendix A and summarized in §3.1, including effects due 
to radiation pressure and anisotropic stress. We show that the location (§3.2) and the heights (§3.3) of the 
peaks as well as all other features in the anisotropy can be simply understood from a physical standpoint in 
this formalism. Moreover the predictions of our analytic treatment are accurate to better than the 5 — 10% 
level in temperature perturbations on all scales, when compared with the full numerical solutions based on 
Sugiyama & Gouda (1992). Finally in Appendix D, we provide a step by step recipe for constructing the 
analytic solutions. 

2. The Evolution of Perturbations 

An exact solution for the evolution of the cosmological perturbations involves solving the coupled evolu- 
tion equations for all of the species present. It is quite apparent that this can only be accomplished through 
numerical integration. However, components such as the neutrinos and the cold dark matter are only cou- 
pled to the photons and baryons gravitationally, whereas before recombination the photons and baryons are 
tightly coupled by Compton scattering. It is thus sufficient to consider the simpler problem of tracking the 
evolution of a single, tightly coupled, baryon-photon fluid in a gravitational potential that accounts for the 
other species. As we shall now show, this problem naturally lends itself to analytic solution for it can be 
described as an oscillator whose restoring force is given by the photon pressure and whose driving term is 
determined by the gravitational potentials. 

2.1 General Equations 

The evolution equation for the fcth Fourier mode of the gauge invariant Newtonian temperature pertur- 
bation* &(ri, ji) is given by the Boltzmann equation with a source from Compton scattering, 

+ ikii{Q + *) = -$ + f [00 - - ^02P2(/^) - il^Vi] , (1) 

where overdots are derivatives with respect to the conformal time rj = J dt(ao/a), fcyU = k -7, with ji as the 
direction cosines of the photon momentum, and f = ^e'^eCTo/ao is the differential optical depth to Thomson 
scattering. Here Xe is the ionization fraction, rig is the total electron density, ax is the Thomson scattering 
cross section, c = 1, and the scale factor a/ao = (1 -|- z)~^. The gauge invariant metric perturbations are \f 
the Newtonian potential, and $ the perturbation to the intrinsic spatial curvature. We will refer to both \f 
and $ as "gravitational potentials." These terms are simply related to the total density fluctuations through 
the generalized Poisson equation and the anisotropic stress (see Appendix A). Note that $ = — \f when 
anisotropic stress is negligible, e.g. in the matter dominated limit. In equation (1), we employ conventions 
where the multipole decomposition is given by 

0(r?,//) = ^(-*f0,(r?)P,(//), (2) 



* For brevity, when discussing a .single k mode of the perturbation, we drop the implicit k index of the 
variables. For example, &(ri, jj,) should be understood as Q{rj, jj,,k). Where confusion may arise, e.g. in the 
discussion of initial and final power spectra, we restore it. Real space fluctuations do not appear in this 
paper. 



and have chosen the amplitude of the baryon velocity 7 • v;, = —ijiV}) accordingly. The appearance of the 
photon quadrupole ©2 in equation (1) merely represents the angular dependence of Compton scattering. 
Finally, we have also assumed a flat geometry. For the open universe generalization of all the arguments 
presented here see Hu & Sugiyama (1994b). A less technical summary of those results may be found in Hu 
(1994). 

Notice that the gravitational potentials \f and $ have two effects on the temperature fluctuations both 
introduced in the original Sachs & Wolfe (1967) paper. The gradient of the Newtonian potential \f induces a 
gravitational redshift on the photons as they travel through the potential well. Since the potential difference 
merely induces a fractional temperature shift of the same magnitude, the combination + \f is the resultant 
temperature perturbation after the photon climbs out of a well of negative \f . We will consequently often 
use + \f to describe the effective perturbation rather than alone. This accounts for what we call the 
ordinary Sachs- Wolfe effect. The time dependence of the metric term $ causes its own time-dilation effect 
referred to here as the integrated Sachs- Wolfe (ISW) effect. 

On the other hand, the baryons evolve under the continuity and Euler equations 

At = -k{Vt-Qi) + ^A^, 

. ^ (3) 

Vi = --Vi, + fc* + f (01 - V^,)IR , 

where R = Zpij/Ap-/ is the scale factor normalized to 3/4 at photon-haryon equality. Here A;, and A^ are the 
baryon and photon energy density perturbations in the total matter rest frame representation (see Appendix 
A). Note that the Newtonian potential acts as a source to the velocity through infall and gives rise to the 
adiabatic growth of perturbations. This coupled photon-baryon system, described by equations (1) and (3), 
fully determines the CMB anisotropics. Notice that the only effect of the other decoupled components is 
through the potentials \f and $. 

2.2 Tight Coupling Limit 

Before recombination, the differential optical depth f is high making Compton scattering extremely 
rapid and effective. Together equations (1) and (3) then imply that Vj, = 0i and &i = 0, for £ > 2. This 
merely reflects the fact that scattering makes the photon distribution isotropic in the electron rest frame. 
Equation (3) then tells us that A;, = 3/4A^, i.e. the evolution is adiahatic. Correspondingly, the resultant 
temperature fluctuations are also adiabatic. This should be distinguished from temperature fluctuations in 
reionized scenarios which are generated by Doppler shifts off electrons at last scattering. In that case, a 
relatively unperturbed photon distribution with Vj, ^ 0i receives a Doppler shift from the last scattering 
event. Here the photons are already isotropic in the electron rest frame, i.e. Vj, = 0i, implying no shift at 
last scattering: the photons merely decouple at recombination. Only where tight coupling breaks down, i.e. 
below the diffusion length, can a Doppler effect arise. However since the term "Doppler peak" is so firmly 
entrenched in the literature, we will continue to use it to describe anisotropics from the acoustic oscillations. 

The tight coupling approximation involves expanding the Boltzmann and Euler equations in the Comp- 
ton scattering time t~^ to eliminate the baryonic variables (Peebles & Yu 1970). To first order, we obtain 
a .single second order differential equation (see Appendix B), 

^o + -j^e>o + k^c^e>o = F(ij), (4) 

a i. -\- II 

where the forcing function F(ri) arises from the gravitational potentials and is given by 



Here the photon-baryon sound speed Cg is 
from which we obtain the sound horizon, 
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Equation (4) tells us that, aside from expansion damping, there are three major tight coupling evolutionary 
effects with different spheres of influence: 

[1] $, the ISW effect on ©o which, when present, dominates at superhorizon scales krj <^ 1; 

[2] k^^, the gravitational infall, which leads to the adiabatic growth of the photon-baryon fluctuations and 
becomes important near the horizon scale fcr? ~ 1; 

[3] k^c^Oo, the photon pressure which cannot be neglected inside the sound horizon kr-g = k J Cgdrj > 1. 

The ordinary Sachs- Wolfe effect, which occurs when the photons stream out of the perturbation after tight 
coupling breaks down, can be taken into account with the combination -|- \f . This partially cancels the 
infall term since it counters the gravitational shift experienced by the photon as it falls into the gravitational 
well. 

The infall and pressure terms must of course compete since the latter prevents the former from causing 
adiabatic growth at sufficiently small scales. The scale at which these forces are in balance is known as the 
Jeans scale. Below this scale, equation (4) describes acoustic oscillations of the photon-baryon fluid. By 
including both the gravitational driving force and the pressure in the intermediate regime crucial for degree 
scale anisotropics, we thus refine the simple Jeans instability argument. 

2.3 Solutions in the Tight Coupling Limit 

Solutions in the tight coupling approximation are straightforward to write down and easily manipulated 
into useful forms (see Appendices B & D). Equation (4) is simply that of an forced, damped oscillator. 
Therefore, the homogeneous F(ri) = equation can be solved by the WKB method, in the limit where the 
frequency is slowly varying. Naturally, the solutions are oscillatory functions with phase (j> = kr-g = k J Cgdrj. 
The WKB approximation is valid for all modes that are in the oscillating regime by last scattering, i.e. 
smaller than the sound horizon at recombination. The particular solution, denoted with an overhat &o, may 
be found from Green's method, 

a/3 1 

[1 + i?(r?)] 1/^00 (»?) = 0o(O) cos fcr,(r?) + ^[0o(O) + -^(O)0o(O)] sin fcr,(r?) 

/3 n ^^"^ 

+ ^ / dr?'[l + i?(r?')]^/^sin[fcr,(r?)-fcr,(r?')]i^(r?')- 
k Jo 

Furthermore, the dipole solution can be obtained from the zeroth moment of equation (1), i.e. the photon 
continuity equation k&i = — 3(0o -|- $). An examination of equation (8) shows that this implies the dipole 
oscillates 7r/2 out of phase with the monopole and has a factor Vg = Cg x (1 + R)~^'^ suppression in 
amplitude. 

The evolution of modes larger than the sound horizon at last scattering can be obtained by taking R = Q 
in equation (8) as proven in Appendix B. The two solutions can be simply joined at kg = 0.08/}^ Mpc~^ if 
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Figure 1. Temperature fluctuation evolution before recombination. The tight coupling 
approximation obtains the evolution to high accuracy compared with the full numerical 
solution once the potentials are known. Here two CDM models are taken as examples. 
Notice that in the high h case, the dipole ©i is significantly suppressed with respect to the 
monopole ©o + ^, and the monopole oscillations themselves are severely modulated. The 
v3 in the amplitude of the dipole accounts for its three degrees of freedom. Here and here 
alone, the numerical results are for a universe which never recombines to eliminate diffusion 
damping at recombination. The arbitrary overall normalization has been set here and in 
the following two figures to \f (0) = —1. 

last scattering occurs sufficiently before photon-baryon equality R(ri^,) ^ 1, as is the case for the standard 
big bang nucleosynthesis scenarios with fi;, ^ fig- 

Already we can gain useful insight on the structure of anisotropics. Scales which reach an extrema in 
the monopole at last scattering will yield a corresponding peak in the anisotropy power spectrum. Note that 
both the positive and negative extrema in temperature yield peaks in the power spectrum. The zeros of the 
monopole will be partially filled in by the dipole but still represent troughs in the final anisotropy pattern. As 
an example, we display in Fig. 1 the tight coupling evolution in the CDM model for two choices of the Hubble 
constant Hq = lOOh km s~^ Mpc~^: (a) h=0.5 and (b) h=0.8. We have also plotted the numerical solution 
to the full set of perturbation equations for comparison. The small discrepancy is in fact almost entirely 



due to slight inaccuracies in the analytic form of the CDM potentials from Appendix A (see also Fig. 3). 
Notice that the suppression of the dipole is more effective in the high h case since R ex Q.},h^ . Furthermore 
the amplitude of the monopole oscillations are severely modulated in this high h case. This feature, due 
to gravitational enhancement of the compressional phase and suppression of the expansion phase, is further 
discussed in §3.3. 

In summary, once the potentials $ and \f are known, the tight coupling solution for the temperature 
perturbation can be readily evaluated with equation (8) at any desired epoch before recombination. In the 
large scale and oscillatory regime, the analytic approximation is limited only by the accuracy with which we 
know the potentials. 

2.^ Diffusion Damping and Recombination 

The first order tight coupling solutions presented in §2.3 can only apply on scales much larger than the 
mean free path of the photons. This scale grows with the expansion and becomes essentially infinite through 
recombination. To handle this photon diffusion problem, we can expand equations (1) and (3) to second 
order in the Compton scattering time t~^. The well known result [see e.g. Peebles (1980) and Appendix B] 
is that the tight coupling solution ©o [equation (8)] is exponentially damped, 

(00 + *) = (00 + *)e-[*/*°('')]% (9) 

where the diffusion scale is essentially the distance a photon can random walk by r], 

fc-2(n) = - / dn ^ , " . (10) 

^ ^'^ 6 Jo 'i (l + i?)2 ^ ^ 

We give the analytic form of t(??), i.e. the ionization fraction Xe(ri), valid through recombination in Appendix 
C. 

This severe damping described by equation (9) arises because as the photons diffuse through the baryons, 
fluctuations become anisotropics which are exponentially damped by scattering (Silk 1968). Here the presence 
of \f represents the ordinary Sachs- Wolfe contribution as mentioned in §2.1. As photons diffuse in and out of 
potential wells, they continue to pick up gravitational redshifts. These contributions are only dependent on 
potential differences and are not damped out by diffusion. Moreover, this treatment automatically accounts 
for the increase in the diffusion length at recombination. As the ionization fraction Xe(ri) —fO,T—fO and 
the diffusion length becomes infinite. Of course, most photons by definition last scatter before the diffusion 
length tends toward infinity. This implies that fluctuations are severely damped under the "thickness" of 
the last scattering surface, i.e. the average diffusion length for a photon at last scattering. Previous work on 
the tight coupling approximation, e.g. J0rgensen et al. (1994), Atrio-Barandela & Doroshkevich (1994), and 
Seljak (1994), have all treated the effects of diffusion and recombination damping separately in a rather ad 
hoc manner. 

The structure of the "Doppler peaks" can be completely described by these damped, driven adiabatic 
oscillations in any scenario where last scattering is sufficiently early. For reionized universes, last scattering 
is delayed, and the diffusion length grows to be nearly the horizon at last scattering. In this case, degree 
scale anisotropics are no longer determined by the adiabatic effect since Vj, ^ ©i but rather by Doppler 
shifts as photons diffuse across the baryons at last scattering. Analytic techniques for studying this situation 
are well known and also accurate at the 10% level in temperature fluctuations (see e.g. Hu, Scott, & Silk 
1994; Hu & Sugiyama 1994a, Dodelson & Jubas 1994 and references therein). These diffusive techniques 
should not be used to describe standard recombination scenarios. 



2.5 The Free Streaming Solution 

Aside from the diffusion modification contained in equation (9), the adiabatic fluctuations present at 
recombination are merely frozen in and free stream to the present (Bond & Efstathiou 1987). They will 
consequently be observable as anisotropics in the microwave background sky today. Let us formalize this 
statement. Equation (1) has the solution (Hu & Sugiyama 1994a) 

[0 + *](??o,y")= /" °{[0o + *-«>V6]f-$ + *}e-^('''''°)e'*''(''-''°)dr?, (11) 

where the optical depth is measured from t] to the present epoch ijo, t(t]i,t]2) = / ^ Tdr], and we have 
dropped the quadrupole term since it vanishes in the tight coupling limit. The combination re"'" is called 
the conformal time visibility function and is the probability that a photon last scattered within dr] of r]. 
Naturally it has a sharp peak at the last scattering epoch r?*- For improvements on the Jones & Wyse (1985) 
fitting formulae for this epoch and recombination in general, see Appendix C. 
Taking the multipole moments and setting Vj, = ©i, we find for £ > 2, 

QiM ~ [00 + *](»?*)(2i' + l)jiikAq,) 

+ QiM[iji-iikAT]^) -{£+ l)ji+i{kAT]^)] .^2) 

+ {21 + 1) / [* - $]jXfcAr?)dr? , 

where Ar? = r/o — r], Ar?* = flo — f]*- The fluctuations on the last scattering surface are determined from the 
aradamped WKB solution (8) as [©o + *](»?*) = [©o + *](??H.)^(fc) and ©i(??h.) = Qi{v*)T^{k), where the k 
index of the perturbations is again suppressed, and we have employed equation (9) to obtain the average 
damping factor, 

Vik)= / °fe-^('''''°)e-[*/*°(''M'dr?. (13) 

Here we have also assumed that all functions save the damping factor are slowly varying compared to the 
visibility function. This is a good approximation on all scales of interest but breaks down for extremely small 
scales where both the oscillations in ji and ©o are rapid. If the potentials are not exactly constant after 
r?H., as is the case if matter-radiation equality occurs close to or after r?*, the integral in equation (12) yields 
an ISW effect after last scattering. Notice that the potential today at the observer does not contribute to 
anisotropics so that the ordinary Sachs- Wolfe effect is given entirely by © -|- \f at last scattering. We do not 
include the present dipole since it cannot be separated from the peculiar motion of the observer. 

Since jt(x) has a peak at £ ~ a;, equation (12) merely represents the free streaming conversion of a 
perturbation on a spatial scale at last scattering to an angular scale on the sky today. An example of 
these last scattering surface fluctuations, we show the analytic results for the CDM scenario in Fig. 2. It is 
important to note that the phase information between the monopole and dipole displayed here is preserved 
in the free streaming transformation. It is not sufficient to free stream the rms sum of the monopole and 
dipole (Bond & Efstathiou 1987) if we want to obtain the detailed structure of the peaks. In particular, 
equation (12) tells us that power in the dipole is distributed more broadly in £ than the monopole. 

Integrating over all k modes of the perturbation, we obtain 

2£+^ V /■dfc fc3|©,(r?o,fc)p 

~1^^'-^JT 2£+l ' ^^^> 

where we have restored the k index implicit in equation (12), and Ci is normalized in the standard manner 
such that the observed (AT/T)^^^ = '^i(2£ + l)CiWi/4:'iT with Wi as the experimental window function. 
We now have the full analytic apparatus to calculate CMB anisotropics in any given model where the 
gravitational potential and recombination history are specified. 
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Figure 2. Analytic temperature fluctuation spectrum at recombination. Again a Harrison- 
Zel'dovich CDM model is chosen as an example. Fluctuations on the last scattering surface 
free stream to the observer creating anisotropics. The phase relation between the monopole 
and dipole as well as their relative amplitudes give rise to the structure of the Doppler peaks. 
Notice that the dipole is significantly smaller than the monopole as expected but is not 
negligible, especially near the zeros of the monopole oscillations. In particular, along with 
the ISW effect, it fills in fluctuations before the first Doppler peak. At intermediate scales, 
gravity is able to shift the equilibrium position of the fluctuations leading to a modulation 
of the monopole peaks (see §3.2). We have drawn in the zero level of the oscillations to 
guide the eye. The kink at fc = 0.04 Mpc~^ is due to the joining of the large and small scale 
solutions. At large scales, the Sachs- Wolfe effect dominates, bringing the effective anisotropy 
to [00 + *](»?*)- V3*(»?*)- 

3. Anisotropies in the CDM Model 

The formalism developed in §2 is applicable to any set of initial conditions or decoupled dark matter 
and may even be readily generalized to cosmological constant or open models. In fact, the open isocurvature 
case is treated in Hu & Sugiyama (1994b). We shall show in this section that it is a powerful and accurate 
technique for calculating and understanding CMB anisotropies by focusing on the standard CDM model. 
This scenario has decoupled cold dark matter supplying the dark mass to make fig = 1 as well as the usual 
photon, baryon and massless neutrino components. In addition, inflation predicts the initial spectrum of 
total density fluctuations to be nearly Harrison-Zel'dovich and adiabatic in the initial conditions (see §A.3 
for the specific definition). For definiteness, we will present the results of this adiabatic Harrison-Zel'dovich 
spectrum. Let us now consider the CMB anisotropies in this CDM model. 

3.1 Gravitational Potentials and Large Scale Anisotropies 

In the spatially flat fig = 1 CDM model we consider here, the potentials $ and \f are constant outside 
the Jeans scale in the radiation dominated epoch and on all scales in the matter dominated epoch. However 
it is important to note that they decay between Jeans crossing in the radiation dominated epoch and full 
matter domination (see Fig. 3). Furthermore, \f = — $ only in the matter dominated limit when anisotropic 
stress becomes negligible. The constant potential approximation, used as a toy model by J0rgensen et al. 
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Figure 3. Gravitational potentials in the Harrison-Zel'dovich CDM model. Notice that 
the potential decays after crossing the Jeans scale in the radiation dominated epoch and 
only flattens out well into the matter dominated epoch. Moreover $ :^ — \f early on due 
to anisotropic stress. These two facts have significant consequences for the temperature 
anisotropics. Notice that the analytic approximations of Appendix A trace the numerical 
potentials reasonably well. 

(1994) and Atrio-Barandela & Doroskevich (1994), is not a good description of the CDM scenario. The 
qualitative difl"erence is that in the CDM scenario, the driving term of the oscillations becomes inefl"ective 
once the Jeans scale grows much larger than the perturbation size. This leads to a more prominent first 
"Doppler peak" and the more complex structure of higher peaks in the real spectrum. 

To understand large scale anisotropics, let us first take the simple zeroth order approximation by ne- 
glecting anisotropic stress and taking the universe to be completely matter dominated at last scattering. In 
this case, the analysis of Appendix A tells us that 0o('?*) ~ —^/s^if]*) ~ —^/s^iflo)- The sign accounts for 
the fact that fluctuations are larger deep in a potential well. The magnitude is given by the initial conditions 
and the ISW efl"ect. Climbing out of the potential \f(??H,) after last scattering, the photons are left with a 
final anisotropy given by [0o + ^]('?*) ~ V3^('?o) as is well known. 

However Fig. 3 shows us that anisotropic stress cannot be ignored in the radiation dominated limit since 
\I>(0) :^ —$(0), and the potentials are not precisely constant until well into the matter dominated epoch. 
The former problem is treated in detail in Appendix A by solving perturbatively for the anisotropic stress. 
The time variation of the potential on the other hand leads to an ISW efl"ect after last scattering and must 
be included if better than 10% accuracy in temperature fluctuations on scales up to and including the first 
Doppler peak is required. This is especially important for low h models where matter-radiation equality 
occurs near recombination as mentioned in §2.5. Because it is an integrated efl"ect, and the angle subtended 
by a given scale increases as time progresses, it tends to spread anisotropics among the low i"s making the 
rise to the Doppler peaks more gradual than it would otherwise be (see Appendix A. 4). This efl"ect and the 
presence of the dipole peak before the first Doppler peak (see §3.2) explains why anisotropics do not behave 
as £(£ + l)Ci = constant for Harrison-Zel'dovich initial conditions as the ordinary Sachs- Wolfe effect would 
imply (Abbott k Schaefer 1986). 
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Figure 4. Individual contributions to the anisotropy in the Harrison-Zel'dovich CDM 
model. At the largest scales (£ < 30), the monopole |0o + ^| from the ordinary Sachs- Wolfe 
effect dominates. The 20% correction from the posi-recombination ISW effect on scales 
larger than the first Doppler peak appears misleadingly small in power (see text). The 
ordinary Sachs- Wolfe effect is overpowered by adiabatic growth of the monopole at small 
scales leading to a deficit at intermediate scales (£ ~ 70) which is filled in by the adiabatic 
dipole ©1 and the ISW effect. Although the dipole cannot be neglected, the monopole 
is clearly responsible for the general structure of the Doppler peaks. Diffusion damping 
significantly reduces fluctuations beyond the first Doppler peak and cuts off the anisotropics 
at £ ~ 1000. 

3.2 Mapping the Anisotropy Spectrum, 

Now let us consider the general features of the anisotropics and determine the scales at which each 
physical process dominates. In Fig. 4, we separate and individually examine the various contributions to 
the anisotropy. In both the large and small scale limits, the monopole ©o + ^ at recombination governs 
the structure of the anisotropics. Scales much larger than the angle subtended by the sound horizon at 
recombination are dominated by the ordinary Sachs- Wolfe effect. However, the correction from the post- 
recombination ISW effect represented in Fig. 4 appears misleadingly small in power. Note that the 20% shift 
in power spectrum normalization from the monopole-only solution is entirely due to the 1% ISW effect. This 
is explained by the fact that the ISW effect adds nearly coherently with the monopole whereas the dipole 
roughly adds in quadrature. Because most of the contribution to the ISW effect occurs near matter-radiation 
equality where r/o — r] ^ r/o, crudely speaking the ISW integral in equation (12) can be estimated as 



[* - $]i,(fcAr?)dr? « [* - $]|;°i,(fcr?o). 



(15) 



The ISW effect therefore mimics the ordinary Sachs- Wolfe effect. In fact, the it partially cancels the ordinary 
Sachs- Wolfe effect at large scales since the potential wells become shallower as time progresses. Indeed, this 
effect contributes at the 20% level in power on all scales up to and including the first Doppler peak for low 
h models. The extent to which equation (15) is a good approximation is discussed in Appendix A. 4. 
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Figure 5. Location of the first Doppler peak. Tlie Doppler peaks are determined by the 
extrema of the monopole oscillations at last scattering h^ = rrn: /rs{rj^,). Because the sound 
horizon becomes independent of fi^ as Q.},h^ goes to zero, the location of the peak is nearly 
independent of fi^ if fi^ ^ fig- The integral over conformal time in r^ makes km oc mh so 
that the angular location of the peaks £„ Ri km^lo ~ mkirjo is nearly independent of h as 
well. 

Near the horizon at recombination, the adiabatic growth of the photon energy density due to infall 
into the potential wells is more than sufficient to make up for the redshift from climbing out. Under the 
sound horizon, photon pressure makes these adiabatic perturbations oscillate and leads to the structure of 
the Doppler peaks. 

The full analytic solution given in equation (8) recovers the location of the peaks to excellent accuracy. 
However, it may be useful to break this solution down to show the intuitive origin of this result. As an 
examination of equation (8) shows, any solution to the tight coupling equations is of the form 



0o(»?) = Ai{vi) cos kr s{vi) + A2{r])s\nkrs{vi). 



(16) 



For adiabatic fluctuations, gravitational infall contributes primarily to the cos kvg harmonic. This is because 
the potential is constant until Jeans scale crossing at which point it begins to decay. The driving term 
therefore mimics cos kvg and causes the monopole to follow suit. This also implies that the dipole goes 
as sinfer^ and vanishes at scales much larger than the sound horizon kvg ^ 1. In contrast, isocurvature 
fluctuations, which begin with vanishing potential, have monopoles which are dominated by the sin kvg 
mode as one would expect. These tendencies hold exactly in the small scale limit (Hu & Sugiyama 1994b). 
Yet even for the first few peaks, it is a reasonable approximation to take cos kvg as the dominant adiabatic 
solution. 

Therefore the peaks in the temperature power spectrum will be located at the scale km which satisfy 
kmfs(ri*) = niTT where m is an integer > 1. Note that this locates the first Doppler peak at roughly the sound 
horizon which is close to, but conceptually distinct from the Jeans scale. The troughs at kf„-i/2^s(f]*) = 
(m — 1/2)17 will be partially filled in by the dipole, including kip which occurs before the first Doppler peak 
(Stompor 1994, see also Fig. 4). In Fig. 5, we plot the scale fci corresponding to the location of the first 
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Doppler peak as a function of ilj, for various choices of h. Since krs{ri) = k J Cgr] = (fc/v3) /(I + R)~^'^dT], 
its dependence on i^i,h^ is weak in the CDM model where fig ^ ^b and i?(??H,) ^ 1. Furthermore, Cg becomes 
independent of h implying that r^ ex r]. Therefore the peak modes are km oc mh. 

Now let us consider the £ space structure of the peaks. Since the last scattering surface is located 
at approximately the horizon distance r/o, the angle subtended by the scale km will roughly correspond to 
f-m ~ kmTio Ri mkirio, and the oscillatory structure of the fluctuations will be partially preserved. Notice 
that the peaks in the dipole are more washed out than the monopole due to the broader nature of the k to 
£ conversion in equation (12) for the dipole. Furthermore, since rj ex h~^, the h dependence of the £ space 
peaks £rn cancels. For the low value of fi;,, which is required by nucleosynthesis, the location of the peaks 
will consequently be roughly independent of both fi;, and h as is well known. Note however that due to 
gravitational suppression of the even peaks (see §3.3) often only the odd m peaks are distinguishable. 

The last significant feature in the spectrum is the diffusion cut off. In Appendix C, we show that the 
optical depth r near recombination is nearly independent of h and only weakly dependent on fi;,. Therefore 
from equation (10), the damping scale becomes ko(ri*) ex /i at recombination. Since the location of the 
peaks scale as km ex h, the cutoff in £ is essentially independent of h and very mildly dependent on fi;,. The 
effect of this damping is displayed in Fig. 4. Employing the undamped solutions ©o and ©i causes a gross 
overestimate of the small scale anisotropics. The problem is compounded in £ space since each £ mode in 
reality gets contributions from a range of k modes, including k > ko- See Appendix C for further discussion 
including the separation of the damping at recombination from that which occurs before it. 

3.3 Heights of the Peaks 

To understand the heights of the Doppler peaks, a somewhat more intricate argument is necessary. Let 
us begin with the first Doppler peak. There are essentially two effects that govern its behavior. Firstly, 
the gravitational infall ex \f competes with the restoring force from the pressure ex cJOq- This would 
naively imply that the deeper the potential, the larger the resultant fluctuation. As we show in Appendix 
A. 4, increasing h makes matter-radiation equality earlier and the potential at the peak deeper due to less 
pressure damping. However anisotropics do not increase with h for all fi;,. A deep potential implies a large 
redshift as the photons climb out after last scattering by the ordinary Sachs- Wolfe effect. Moreover if the 
sound horizon is approximately the particle horizon r-g = J Cgdrj ^ rj, the window in which adiabatic growth 
dominates is small (see §2.2). Thus, a deeper potential implies a sm«//er final anisotropy if the pressure is 
held fixed at a high value, i.e. Cg ~ l/v3- 

Therefore the way to increase the anisotropics at the first peak is to reduce the pressure, i.e. the sound 
speed, rather than deepen the potential. This can be accomplished if we raise Q.},h^ . The height of the first 
Doppler peak thus increases with fi;, (Holtzman 1989, Fukugita, Sugiyama, Umemura 1990). As for the h 
dependence, we have uncovered two opposing effects. Since the sound speed becomes Cg ~ l/vS, independent 
of h as fift goes to zero, the variation of the pressure with h asymptotically disappears. This is exactly the 
limit where the Sachs- Wolfe mechanism is most effective at counterbalancing adiabatic growth. Therefore at 
low fift raising h will decrease the anisotropy whereas at high fi;, it will increase the anisotropy. The crossover 
point happens to coincide approximately with the nucleosynthesis value of fi;, Ri 0.05 if /i Ri 0.5 — 0.8. Thus 
the h dependence of the first Doppler peak for standard CDM will be relatively weak. Note however that the 
Sachs- Wolfe mechanism becomes less important at smaller scales due to the decay of the potential, whereas 
the pressure argument is scale free. 

Now let us consider the general structure of the peaks. For all odd numbered peaks, the situation will 
be qualitatively identical to that of the first Doppler peak: \f and ©o have opposite signs corresponding 
to overdensities in potential wells and underdensities at the peaks. In other words, inside the gravitational 
well, the oscillation is in its compressional phase and is thus enhanced by gravitational infall. On the other 
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Figure 6. The CDM anisotropy spectrum. Notice that for high fi^ raising h increases 
the amplitude of the peaks, due to a lowering of the pressure, whereas for low fi;, raising h 
decreases it, due to a deepening of the potential well out of which the photon must climb. 
The pressure argument also explains the monotonic increase in the heights with fi;,. The 
high fift/i'^ case also has a significant decriment between the first and third Doppler peaks 
due to the modulation eff"ect and a significantly smaller dipole contribution. 

hand, the even numbered peaks correspond to expansion inside the well and are correspondingly suppressed. 
The more gravity dominates over pressure, i.e. the higher fi;,/}'^ is, the more eff"ective is this pattern of 
alternating enhancements and suppressions. In some cases, the second Doppler peak disappears altogether. 
This modulation eff"ect can also be thought of as a shift in the equilibrium point of the oscillations due to 
the gravitational force (c./. Fig. 2). Thus it disappears at high h where the potential becomes negligible. 
Since in the high fi;,/}'^ case the dipole is also severely suppressed with respect to the monopole peaks, this 
implies that in these models the decrement between the first few positive (odd) peaks will be substantial. 

In Fig. 6, we display the final anisotropy spectrum Ci for various choices of fi;, and h. Notice that for 
fij = 0.06, raising h increases the amplitudes of the peaks, whereas for fi;, = 0.03 it decreases the amplitudes 
of the peaks. Furthermore, the modulation eff"ect described above shows up for the high ^j^h^ case making 
the even peaks invisible. In all cases, the analytic solution matches the full numerical results typically to 
better than 10% in power. Our analysis therefore includes all the major contributions to the anisotropy. 
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4. Conclusions 

We have presented a fully analytic treatment of CMB anisotropies under the assumption of tight coupling 
between the baryons and photons before recombination. Our technique describes the tightly coupled photon- 
baryon fluid as an oscillator in an arbitrary potential well. Since dark matter affects this system only through 
these potential wells, all dark matter models can be treated under this formalism. Applying this method to 
the CDM model, we show that it typically obtains CMB anisotropies to 10% or better in power on all scales. 
Moreover it allows us to separate and explain each and every process that contributes to the final anisotropy. 
On the very large scale, the familiar Sachs- Wolfe tail dominates the CMB spectrum. However the flatness of 
this tail is broken by both the adiabatic photon dipole contribution and the ISW term, particularly for small h 
models because the universe is not totally matter dominated at last scattering surface. The series of "Doppler 
peaks" results from extrema in the oscillatory monopole at last scattering. These acoustic oscillations are 
driven by the gravitational potentials which cause the adiahatic growth of density fluctuations. In fact the 
name "Doppler peaks" is itself somewhat misleading since it implies that they result from Doppler shifts 
induced by last scattering off electrons in infall. Since anisotropies in a reionized scenario are generated by 
such a mechanism, it is perhaps preferable to call these oscillatory fluctuations "adiabatic peaks" instead. 

We moreover obtain the location and heights of these peaks to excellent accuracy and show that they 
are easily explained by variations in the sound speed and the depth of the potential wells. Thus we have 
also clarified the dependence of CMB anisotropies on cosmological parameters, i.e. the Hubble constant h 
and the baryon fraction fi;,. Furthermore, by taking the recombination process into account more properly 
than previous work (see Appendix C), we recover the damping scale from photon diffusion arguments alone. 
This implies that the damping due to the finite thickness of last scattering surface and photon diffusion are 
the one and the same. 

Our interpretation of these physical processes is confirmed by its excellent agreement with the full 
numerical solutions of the coupled evolution equations. Readers can even use our method, reconstructed 
in detail in Appendix D, instead of full Boltzmann code if only 10% accuracy in power is required. In 
fact, our formalism is useful in combination with numerical work as well. It can be employed as a simple, 
portable standard with which various numerical treatments can be compared. Moreover, in the large scale 
and oscillatory limit, the approximations presented here are limited only to the extent that the gravitational 
potentials are known. Since the evolution of density fluctuations can be solved far more easily than the 
infinite heirarchy of temperature multipoles, our formalism can be used as a quick and simple way to 
accurately calculate and understand the CMB anisotropy in any model where the fluctuations are still linear 
at recombination. 
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Appendix A: Gravitational Potentials and Anisotropic Stress 

A.l General Relations 

The potential $ is related to the total density fluctuation Ay = Spx/pT by the generalized Poisson 
equation 

fc2$ = 4TrGp{a/aofAT, (A-1) 

and thus accounts for the decoupled cold dark matter and neutrinos. Throughout this appendix, we normalize 
the scale factor a at matter-radiation equality. In order to determine the potential, we must analyze the 
contributions from all species. In fact, the evolution equations for the CDM and the massless neutrinos can 
be obtained from equations (1) and (3) by the replacements A;, —f Ac and & —f N with t —f where Ac and 
N are the density fluctuation in the CDM and the temperature perturbation of the neutrinos respectively. 
These replacements are valid since the only difference between these species in the non-degenerate limit is 
Compton scattering. 

We also need the Newtonian potential \f which is related to $ through the anisotropic stress 11, 

where 11 is given by the quadrupole moments of the the photons and neutrinos, 

pl[ = —(p.,e2 + P.N2). (A-3) 

If either the pressure p is unimportant or the quadrupole moments vanish, \f = — $. We know that ©2 
vanishes due to the isotropizing effect of Compton scattering in the tight coupling limit. Even N2 can only 
be generated through free streaming from Nq. Our approach to obtaining an analytic solution for $ and \f 
is therefore to include 11 as a small perturbation. 

Combining the four evolution equations, we obtain (Kodama & Sasaki 1984) 

At - 3w-At = -(1 + w)kVT - 2-wIl, 

"4 "9 (A-4) 

Vt + -Vt = -TrkAr -I- fc* k 11, 

a 3(1-Fw)2 ' 3 1-Fw ' 

where w = p/ p, and the total velocity is given by the sum over component velocities (p+p)Vt = '^f(pi+Pi)Vi, 
with Ky = ©1 and V^ = Ni. Here we have assumed adiabatic conditions and have also used the relations 



A^ =400 + 4 



a Vt 
a k 



A, = 47Vo + 4--^, 
a k 



in converting from the Newtonian variables to the total matter rest frame variables. 
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A. 2 Zeroth Order Large Scale Solution 

In the zeroth order approximation where H = 0, equation (A-4) has exact solutions in the limit that 
we can ignore the fcAy pressure term in the velocity equation. This is appropriate for superhorizon sized 
fluctuations.* The growing mode AT(a) = AUaia) is given by (Kodama & Sasaki 1986), 



Uaia) 



a- + -a^ 



—a 
9 



16 16 , 

1 Va + 1 

9 9 



a(a + l) ' 



(A-6) 



where recall that the scale factor is normalized at matter-radiation equality Ugq = 1. Here the normalization 
factor A{k) defines the relative weighting of the k modes through the initial power spectrum (see A. 3). On 
the other hand for reference, the decaying mode takes the form 



Uoia) 



a\/a + 1 



(A-7) 



Note that Uq = (10/9)a'^ in the radiation dominated limit and a in the matter dominated epoch. This gives 
the well known result that outside the horizon, the potential is constant in both the radiation and matter 
dominated limits. In this 11 = limit, we obtain \f (??o) = ^/w^{Q). Equation (A-6) with (A-4) and (A-5) 
also implies that 0o(O) = - V2*(0). 

This behavior for the potentials of course has significant consequences for CMB anisotropics. Equation 
(1) tells us that at large scales where the ISW effect dominates, ©o = — ^- This implies that the fluctuations 
at last scattering, which determine the ordinary Sachs- Wolfe effect, are given by 



0o(»?*) = 0o(O) + *(r?*)-*(O), 



(A-8) 



since $ = — ^. If last scattering occurs well into the matter dominated epoch, the potential is constant and 
^(r^H.) = ^{rjo). Putting these results together, we obtain 0o('?*) = — V3^('?*) ^^^ ^^e total Sachs- Wolfe 
temperature perturbation after the photons climb out of the potential. 



[0o + *](r?*) = ^*(r?o). 



(A-9) 



which is the familiar Sachs & Wolfe (1967) result. Indeed, this zeroth order solution can in fact be employed 
at large scales if only 10% accuracy in the temperature fluctuations even at the COBE DMR normalization 
scale is acceptable. Of course, an error at the normalization scale causes an error on all scales. It is therefore 
preferable to correct for the small effect of anisotropic stress perturbatively. 



* Since inside the Jeans scale, pressure acts to damp the growth of perturbations, we will denote the 
■undamped large scale solution with an overbar analogous to our notation for the tight coupling solutions. 
This is to caution the reader that the solutions only apply prior to horizon crossing for all modes that cross 
before matter-radiation equality. 
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A. 3 First Order Large Scale Solution 

Before recombination, due to the isotropizing effects of scattering, the anisotropic stress of the photons 
is negligible small. Hence the main contribution to H comes from the neutrino quadrupole anisotropy N2 
[see (A-3)]. Here we analytically obtain the growing mode solution for density perturbations including the 
contribution of anisotropic stress for modes outside the horizon at matter-radiation equality. This together 
with the matter transfer function is sufficient to obtain the gravitational potential on all scales. 

We take into account the neutrino quadrupole anisotropy perturbatively. Namely, we use the exact 
zeroth order solutions (A-6) and (A-7) to obtain the anisotropic stress. We then take this solution to 
iteratively correct for anisotropic stress in equation (A-4). If we neglect higher order multipole components, 
which is reasonable for superhorizon sized modes, the second moment of the the Boltzmann equation for the 
neutrino becomes 

N2 = hNi « hvr , (A-10) 

where we employ causality to infer that the fluid cannot separate at superhorizon scales, i.e. Ni = V^ ^ Vp 
(see Hu & Sugiyama 1994b). Moreover the exact zeroth order solution for Vp is obtained using equation (A-6) 
in the first of equations (A-4), yielding the solution to equation (A-10), 

- / X r da' 1 f , , , ,dVG\ 

NM - 2a/ ^^7^ \^o - («' + 1)«'^) , (A-11) 

where we have used the relation Zw = 1/(1 + a) and recall that the overbar represents the aradamped 
superhorizon solution. Although it is possible to analytically integrate equation (A-11), the expression is 
cumbersome. Instead, we can employ an approximate solution which is exact in the limit a <^\ and a ^ 1, 

- , , 1 20a + 19 8 a 8, /3a + 4\ 

N2(a) = — AVg A+-ln — ^\A. (A-12) 

^ ' 10 3a + 4 3 3a + 49V4y ^ ' 

We have checked that this approximation works quite well by comparing it to equation (A-11) and the full 
numerical solution. 

Next, we employ the above solution for N2 in equations (A-4). These two first order equations may be 
rewritten as one second order equation for Ay. Analogously to the tight coupling equation (4), the particular 
solution including the source terms H and H can be obtained from the homogeneous solutions Vq and Vd 
by Green's method, 

AT(a) = f 1 + '^-f^ AVaia) + ^/.[Ii(a)?7G(a) + h{a)Vo{a)\ , (A-13) 



where Ii(a) = Q da' S{a')Vo{a'), hia) = J^ da'S(a')UG(a'), 



24 (a + 1)5/2 r 2a d 
2 



5(a) =- ^"7^^- /^^^A^^(a) 



5 a2(3a + 4) [ 3a + 4 da 

AUaia) 



(3a + 4)(a+l)' 



12 



(a +1)2 a+1 3a- 



(A-14) 
N2(a)j, 



and /jy is the ratio of neutrino to total radiation density f^ = Pv/ifiv -\- P~i)- If we assume three massless 
neutrinos and the standard thermal history, Pv j p~f = 3(7/4)(4/ll)'*'^/2 = 0.68, i.e. f^, = 0.405. The first 
term in equation (A-13) comes from the initial conditions for Ay which can be iteratively established by 
employing equation (A-10) in (A-4). All terms which are proportional to f^ in the equation (A-13) come 
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from equation (A-3) since the anisotropic stress H Ri (12/5) f^N2. The asymptotic behavior of the equation 
(A-13) is 



At (a) 



l + -U]AUG(a) (a<l) 



(A-15) 



1 + -Ml- 0.333) jAUaia) (a > 1) . 



Here we have used the fact that if a ^ 1, the decaying term I2U0 may be ignored and /i —f —0.333. 
Therefore we may obtain a simple approximate expression for the large scale density fluctuations, 



At (a) 



2 . / a 

1 + -/J 1-0.333^ 
5 V a+ 1 



AUG(a). 



(A-16) 



Again we have checked this approximation works reasonably well by comparing it to numerical calculations. 
The potentials $ and \f are therefore written as 



#(a) 



3 fk. 

4 



-AT(a), 



«(«) = -I (^ 



1 /. . .^8 N^ja) 
AT(a) + -u — — r 



(A-17) 



where k^q = v2(f^o-ffoao)^ is the scale that passes the horizon at matter-radiation equality. By using the 
asymptotic form of At and N2, we easily obtain the corresponding relation between $ and ^, 



#(a) 



! + -/.)*(«) (a<l) 



(A-18) 



-*(a) (a > 1). 



Also of interest are the ratios of initial to final values of the the gravitational potentials: $(ao) = — ^(ag) = 
0.86$(0) and \f (ag) = 1.00\f (0). Thus we see that the correction for anisotropic stress makes a 10% difference 
in \f in the radiation dominated epoch. If recombination occurs near equality, this results in a small correction 
to the standard Sachs- Wolfe formula due to anisotropic stress. 

The initial conditions for the perturbations may now be expressed in terms of $(0), 



*(0) = *(0) = -0.86$(0), 
0(0) = 0(0) = 0.43#(0). 



(A-19) 



Note that since all modes are superhorizon sized at the initial epoch, the overbar is superfluous. Moreover, 
even in the initial conditions, the anisotropic stress represents a small but important correction to the H = 
solutions of §A.2, $(0) = — \f (0) = 20o(O). Finally, we can relate these quantities to the initial power 
spectrum. 



fc3|$(0,fc)p = fc3|#(0,fc)|^ 



5 / 2 , 



-^ ) A'ik) = Bk"-^ 



(A-20) 



where we have again restored the implicit k index and the overall normalization factor B is fixed by the 
COBE DMR detection (see Appendix D). The Harrison-Zel'dovich initial spectrum predicted by inflation is 
obtained by setting n = 1. 
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A. 4 Small Scale Potentials 

Next we need to obtain solutions of \f and $ in the small scale limit where pressure cannot be neglected. 
Qualitatively speaking, we know that the potentials decay inside the Jeans length in the radiation dominated 
epoch since pressure prevents Ay from growing. However in general, it is impossible to obtain the exact 
solution valid through matter-radiation equality even if we neglect the anisotropic stress term. Only the 
asymptotic behaviors in certain limits have been found (Kodama & Sasaki 1987). For the CDM scenario, it 
is well known that the final value of the potential at small scales is obtained from the superhorizon solution 
(A-17) by the transfer function $(ao) = -*(ao) = T(fc)#(ao), where 

Tik) = Ml_L^^[i + 3.89Q + (UAqf + i5A6qf + i6.71q)T"' , (A-21) 

with q = fc/[fio'j'^ exp(— 20;,)] (Peacock & Dodds 1994, Bardeen et al. 1986). Note that q ex k/keq ap- 
proximately, reflecting the fact that only modes that cross the Jeans length before equality are suppressed. 
This implies that the potentials are larger in amplitude if equality occurs later, i.e. for high i^oh^ models. 
Equation (A-21) therefore empirically accounts for the lack of growth in the radiation dominated era. 

Now let us consider the time evolution of the potential. We know that in the matter dominated epoch the 
potentials are constant on all scales. Therefore, we smoothly join the superhorizon scale solutions of equation 
(A-17) with a constant matter dominated tail whose relative amplitude is given by the transfer function. 
Since the Jeans crossing epoch is approximately same as horizon crossing time in radiation dominated era, 
we can take {k/ Ha) ~ ak/keq ~ 1 as the matching epoch, 

$(a) =#(a) |[1 - Tik)] exp[-ai (ak/k.^)^] + T(fc)| , 

(A-22) 
*(a) =*(a) {[1 - T(k)] exp[-a2 (ak/k.^)"] + T(fc)} , 

where ai, a2 and /3 are fitting parameters. We also need a small correction to take into account the free 
streaming oscillations of the neutrino quadrupole inside the Jeans scale. A very simple approximation can 
be obtained by making the replacement ^2(0) —^ N2(a)cos[0.5k/(Ha)] in equation (A-17) for \f(a). Here 
0.5 is a fitting factor and the Hubble parameter H(a) = (a/a)(ao/a). This crude approximation is sufficient 
for our purposes. Comparing this functional form (A-22) with numerical results, we obtain a good fit for 
ai = 0.11, a2 = 0.097 and /3 = 1.6. 

In order to calculate the post-recombination ISW effect, we take the direct derivative of equations (A-22). 
Although this makes the estimation of the derivatives much worse than for the potentials themselves, it is 
a reasonable approximation in the context of temperature fluctuations since the ISW effect is effectively 
only a perturbation to the spectrum. Another way of seeing this is to compare the full calculation to the 
approximation to the ISW integral, 

r [* - ^]u(kAv)dv «[*-$] HlJiikvo). (A-23) 

■Jv 

This simple approximation eliminates the need for detailed knowledge of the derivatives by assuming that 
most the contribution comes from early enough epochs that ji(kAri) Ri ji(krio). Since there is only a 
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Figure 7. The ISW effect. Ignoring the ISW effect entirely leads to a significant error 
in both the normalization at 10° and shape of the anisotropics due to contributions near 
recombination. This can be partially accounted for by approximating all of the ISW con- 
tribution to occur near recombination. This approximation leads to 10 — 15% errors in 
temperature due to the fact that some of the contribution comes from more recent times 
where the fluctuation subtends a larger angle angle on the sky. The full integration therefore 
has more power at larger angular scales and makes the rise to the first Doppler peak more 
gradual. 

contribution near matter-radiation equality, the integral does indeed get most its contribution at early 
times. 

In Fig. 7, we show the results of ignoring the ISW integral, using equation (A-23), and integrating 
the ISW effect with the derivative of equation (A-22). Ignoring the ISW contribution entirely is clearly 
a bad approximation at the 20 — 25% level in temperature fluctuations since it seriously misestimates the 
normalization and the shape of the rise toward the first Doppler peak. Employing the approximation improves 
the situation to the 10 — 15% level in temperature. This approximation has the effect of putting the power 
at too high a multipole £ corresponding to the angle that the scale subtends at last scattering rather than 
the true distance. We see that detailed knowledge of the form of the derivatives only represents a 10 — 15% 
shift in power which justifies our crude approximation of it. 

Appendix B: The Tight CoupUng Approximation 

B.l The WKB Approximation 

The tight coupling approximation has often been used in the past to describe the behavior of fluctuations 
at sufficiently small scales that gravitational effects can be neglected entirely (e.g. Peebles & Yu 1970). Here 
we show that it is easy to include the effects of a realistic time dependent potential in the tight coupling 
formalism even in the intermediate regime where both gravity and pressure play a role. In short, the 
technique involves expanding the Boltzmann equation (1) and Euler equation (3) in the Compton scattering 
time T~^ in the limit that Compton scattering is much more rapid than expansion or gravitational infall. 
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To zeroth order, we regain the tight coupling identities, 

0o(»?) = 0(»?,/^) + «>0i(»?), 

where we have used the fth moment of equation (1) to show that 0£ = for £ > 2. These equations 
merely express the fact that the radiation is isotropic in the baryon rest frame. Substituting the zeroth order 
solutions back into equations (1) and (3), we obtain the iterative first order solution, 

0o = --0i-*, 

• (B-2) 

1 + R 1 + R ' 

where we have used the relation R = (a/a)R. Note that we have used the tight coupling approximation to 
eliminate the multiple time scales and the infinite hierarchy of coupled equations of the full problem. In fact, 
this simple set of equations can readily be solved numerically [see e.g. Seljak (1994)]. This may be desirable 
in the case where one needs to know very accurately the transition regime between the large and small scale 
analytic approximations. 

For the purpose of obtaining analytic solutions, it is preferable to rewrite equation (B-2) as a single 
second order equation, 

00 + 7^00 + k^c^Oo = F(ij), (B-3) 

where 

F{V) = -* - Y^* - y*, (B-4) 

is the forcing function with $ as the ISW effect, $ as the modification to expansion damping, and \f as the 
gravitational infall. The homogeneous F(ri) = equation yields the two fundamental solutions under the 
WKB approximation, 

0Jt]) = (1+ m-i/^cosfcr,, 

0iiq) = il + R)-^'\mkr,, 

where the sound horizon r^ = J Cgdr] can be integrated to give 



3 V R[Veq) 1 + ^/R(Veq) 

with kgq = v2(f^o-ffoao/ae5)^ — {^ ~ '^y^)l'f]eq as the scale that enters the horizon at equality. The phase 
relation just reflects the nature of acoustic oscillations. If the sound speed were constant, it would yield the 
expected dispersion relation w = kcg. 

Now we need to take into account the forcing function F{r]) due to the gravitational potentials \f and 
$. Employing the Green's method, we construct the particular solution, 

6„(„ = CM,) + CM,) + /' ;;";'^;";!~-i''!^';'V c'w- (b-^) 

■h Oa(r]')0b(r]') - 0a(r])^b(r]) 
Employing equation (B-5) yields 

Oaiv'Mv) - OaivMv') = [1 + R{V)]-^'^[1 + R{v')]-^'^sm[krs{v) " kr,{v')] , (B-8) 

and 

OaWWiW) - 0aW)0iW) = ^[1 + RWT' ■ (B-9) 

Equation (8) now follows by fixing Ci and C2 with the initial conditions. 
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B.2 Large and Small Scale Corrections to the WKB Approximation 

On large scales the WKB approximation breaks down, whereas on small scales we need to modify 
equation (B-2) to account for photon diffusion. The WKB approximation assumes that the frequency, i.e. 
the sound speed, is not rapidly varying 

ikc,f >(1 + R)"''^ii + R)-^l\ (B-io) 

and is valid if i? ^ 1 and on small scales. Roughly speaking, this requires the mode to be in the oscillatory 
region by last scattering. More specifically, it is reasonably well satisfied at recombination if kg > O.OSh^ for 
the range of fi;, consistent with nucleosynthesis. Since we know the solution at scales k <^ kg where we can 
neglect pressure, we can obtain the full solution by matching the two. Let us see how this is done. 

At large scales, we know the only the ISW effect is important since f-mode coupling and infall are only 
effective inside the horizon. The solution to equation (B-2) is therefore 

0o(r?) = 0o(O) - / $(r?')dr?'. (B-11) 



Jo 
On the other hand, the WKB approximation at large scales, i.e. the krg —f limit, predicts 

(1 + Rf^Qoi^) = 0o(O) - y H^')^^ {[1 + RW)f'^k[rg{^') - rg{^)]}d^', (B-12) 

where we have integrated once by parts and employed 0o(O) = —$(0). If i? = 0, r^ = c^ = l/v3 and the two 

expressions are identical. Therefore, the large scale solution obeys equation (B-7) if the two fundamental 

solutions are taken to be 

0a(v) = cos kr-g. 

(B-13) 
0b(r]) = sinkr-g, 

This is a particularly useful form to express the large scale solution because in the CDM model, R(ri^,) <^ 1 
since fig ^ ^b- This solution will therefore approximately join into the WKB approximation at small scales. 
We take the matching point to be kg = O.OSh^ but the results are not terribly sensitive to what scale is 
chosen (c.f. Fig. 2). In fact for a simple approximate estimate (~ 20% in AT/T) to the CDM spectrum, 
one can employ equations (B-13) on all scales. More explicit formulae can be found in Appendix D. 

A correction to the tight coupling equation (B-3) itself must be applied at small scales. Notice that 
the finite time scale of Thomson scattering, controlling the mean free path of the photons, has not appeared 
anywhere in the analysis so far. That is because we have only expanded to first order in t~^. For diffusion 
effects, we need solve iteratively to second order. In this case, we do neglect the effects of the ISW term 
and gravitational infall for the baryons. The dispersion relation must consequently be modified to have an 
imaginary term. If we write the solution for Oq + \f as exp(i(/)£)) we obtain [see e.g. Peebles (1980)] 

4>Diil) = ±k f cgdq + i[k/kDiil)]\ (B-14) 

Jo 

where the damping length is given by 

, 2/ X 1 /■ , 1 -R^ + 4(1 + i?)/5 



25 



u 

SIX 



I Q 




0.1 - 



0.01 



Figure 8. The evolution of damping through recombination. The damping length increases 
rapidly at recombination corresponding to the fact that the mean free path of the photons 
becomes infinite. Weighting by the visibility function, which tells us when the photon last 
scattered, yields the average damping factor. We have also included a simple approximation 
to the damping integral which can be used for estimation purposes before recombination. 

This then yields the damping correction given in equation (10). Before recombination, this integral may be 
simply approximated as 



fc-2 ^ 1.7 X 10^(1 - Yp/2)-H,^bh^)-H.^oh^)-^'^ (—) 



5/2 



^^Oeq/a + 2 



Mpc^ a<^a^, (B-16) 



where Yp Ri 0.23 is the primordial helium mass fraction. This may be useful for estimation purposes but 
should not be used to describe the detailed damping process at recombination since the damping length 
suddenly increases to infinity. To illustrate this effect, in Fig. 8 we plot the evolution of the damping 
scale through recombination obtained by following the true ionization history obtained in Appendix C in 
comparison with equation (B-16). 

Appendix C: Recombination 

Following Peebles (1968) and Jones & Wyse (1985), we solve for the ionization history through recombi- 
nation. Since we wish to obtain the detailed behavior of photon diffusion damping including its dependence 
on cosmological parameters, we need to improve upon the fitting formulae obtained by Jones & Wyse (1985). 
This is especially necessary in CDM scenarios where fi;, ^ fig- Note that the full numerical treatment em- 
ploys the numerical values for the ionization history rather than the approximations presented here. 

The total optical depth from the present to the critical recombination epoch 800 < z < 1200 can be 
approximated as 

<^^o)^nl^(-^)" (c-i) 
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Figure 9. The redshift visibility function. Notice that the weak dependence on fi;, of the 
visibility function is adequately described by the analytic fitting formula, whereas the Jones 
& Wyse (1985) fitting formula [their equation (23)] does not. 

where ci = 0.43 and C2 = 16 + l.Slnfi;,. Since the range of reasonable values for h is limited to be 
0-5 ^ h < 0.8, we have ignored the small h dependence. For definiteness, we take last scattering to occur at 
Zh, where the optical depth t(z^,, 0) = 1. It immediately follows from (C-1) that this occurs at 
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-Ci/C2 



fi, 



-0. 027/(1+0. llln Hi,) 



which is weakly dependent on fi;,. The differential optical depth f then becomes 



f(z) 



C2 

1000" 



■^l 



1000 



C2-1 



-(1 + z), 



(C-2) 



(C-3) 



where f is by definition positive since f = d[T(ri', ri)\/drj. Finally the ionization fraction is given by Xe{z) = 
Tao/ne(JTa where 

(rieO-Ta/ao)"^ = 4.3 x 10^(1 - Yp/2)-^{nih^)-^{l + z)-^Mpc (C-4) 

with Yp Ri 0.23 as the primordial helium mass fraction. Of course, where the formula (C-3) implies Xe > 1, 
set Xe = 1, i.e. t = rieCTa/ao- In Fig- 9, we show the numerical values for the visibility function in redshift 
space — (dr/dz)e~'" compared with these analytic fits. 

The decrease in ionization fraction implies an increase in the distance a photon can diffuse and damp 
(see Fig. 8). This is evaluated by employing the above formulae for the ionization fraction in equation (10). 
Recombination therefore causes a sudden increase in the damping of anisotropics. In Fig. 10, we display the 
effect of recombination on anisotropics by comparing the full solution for the damping with that obtained by 
assuming instantaneous recombination, i.e. approximating the damping factor as T>(k) k, exp{ — [fc/fc£)(r?H.)]'^} 
where k~jj [t]) is given by equation (B-16). The misestimation of the damping scale is significant, but it is 
still a much better approximation than neglecting damping entirely. 
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Figure 10. The effect of the finite thickness of the last scattering surface. Estimating the 
damping in the instantaneous recombination approximation leads to a significant underes- 
timate of the damping scale. It is however far better than neglecting diffusion damping 
entirely. 

Appendix D: A User's Manual 

D.l Explicit Tight Coupling Expressions 

In this Appendix, we bring together and list the various ingredients necessary for the analytic solution 
that have been scattered throughout the text and other appendices. First of all, although equation (8) is the 
best way to write the solutions in the tight coupling limit for understanding the physics of CMB anisotropics, 
for calculational purposes it is convenient to express the solutions in a more explicit but cumbersome form. 
One advantage of the analytic tight coupling solutions is they do not require the use of time derivatives of 
the potentials despite the appearance of equation (8). Thus accuracy is not compromised by our lack of a 
detailed description for $ and \f . Integrating equation (8) by parts twice, we obtain 



(1 + i?)i/4[0o(»?) + *(»?)] = [cos kr,(r]) + J(0) sin fcr,(r?)] [0o(O) + $(0)] + I(r]) , 



where the overhat denotes the undamped solution 

k drj 



J(,)^_(l + JJ)3/4V3i.(l + JJ)-l/4 a/3 R 



4fc yiT^' 



and 



with 



I{r]) = —= \ dr]' ^{r]')G{r]') sm[krs{r]) — krs{r]')] 
V 3 Jo 



G(v) = (l + R) 



-1/4 



$ 4fc2 



(D-1) 



(D-2) 



(D-3) 



(D-4) 



Here we have employed the identity 0o(O) = —$(0). Since the ISW effect predicts constant ©o + ^ at 
superhorizon scales, we have written these expressions in terms of that quantity. 
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The dipole solution ©i can be similarly obtained from the photon continuity equation k&i = — 3(00+^), 

(1 + Rf''^^ = [1 + J('?)J(O)][0o(O) + $(0)] sinfcr,(r?) 

+ [Jiq) - J(0)] [0o(O) + $(0)] cos fcr, (r?) (D-5) 

+ J{r])I{r]) -^ dr]^{r]')G{r]') cos[krs{r]) - krs{r]')] , 

where we have used the relation Vg = Cg = (l/v3)(l + R)~^'^ . Notice that we do not need $ even in the 
boundary terms in either equation (D-1) and (D-5). These forms also bring out the fact that whereas the 
monopole term is ex (1 + R)~^''^, the dipole being ex (1 + R)~^''^ is suppressed in comparison. 

On the other hand, the large scale solution discussed in Appendix B may be obtained by dropping R. 
To be explicit, these are 

[0o(r?) + $(r?)] = [0o(O) + $(O)]cosfcr,(r?) + ^ / dr?'[$(r?') "*(»?')] sm[fcr,(r?) - fcr,(r?')] > (D-6) 

and 

0i(,;)/V3= [0o(O) + $(O)]sinfcr,(r?)-^ / dr?'[$(r?') "*(»?')] cos[fcr,(r?) - fcr,(r?')]- (D-7) 

V 3 Jo 

Finally, the following relations are useful for computation: 

R=—^-J^a, R = aR{7jeq) = ^VlT^Ri^eq), R=-k%R{7ieq), (D-8) 

i — J^ 4 \Iq \J I 4 

where we have employed the relation k^qf] = 2v2(\/l + « ~ !)• Here 1 + Pv j p~f = (1 — fv)~^ = 1.68. Note 
that the scale factor is normalized at equality Ggq/aQ = a^ = 4.0 x 10~^ (Tq /2.7K)'^(iloh^)~^ , and the scale 
which passes the horizon at equality is keq = l.n/rjeq = 7.46 x 10~^(To/2.7K)~^i^oh^ Mpc~^, with To as 
the present temperature of the CMB. We have gathered together many of these commonly used symbols, 
and the equations in which they are defined or first appear, in Table 1. 

D.2. Poor Man's Boltzniann Code: A Recipe 

Now let us outline the steps in the analytic calculation: 

[1] For k > kg = 0.08/}^, take the tight coupling solutions for the undamped monopole 0o and dipole 0i 
from equations (D-1) and (D-5) with the potentials from equation (A-17) and (A-22) and the initial 
conditions from equation (A-19) and (A-20). Evaluate this at last scattering r?* given by equation (C-2). 

[2] For k < kg = 0.08/}^, repeat the steps in [1] using the large scale solutions equations (D-6) and (D-7) in 
place of equations (D-1) and (D-5). Join the two solutions at kg. 

[3] Evaluate the damping scale function ko(ri) from equation (10). Use the recombination fitting formula 
(C-3) for f where it implies the ionization fraction Xg < 1 and equation (C-4) for earlier epochs, i.e. 
Xe = 1 and f = rieCTa/ao- Integrate the damping factor against the visibility function -rexp(— r) from 
equation (C-1) and (C-3) thus calculating T>{k) from equation (13). This damping factor is defined such 
that the true anisotropy at last scattering [0o + \f](??H,) = V{k)[Qo + \f](r?H.) and 0i(??h,) = V{k)Qi{rj^,). 

[4] Free stream the perturbation from last scattering to the present by employing equation (12). For an 
added 10% correction in the low h CDM scenarios, perform the ISW integral in equation (12) using the 
derivative of the potentials in equation (A-22). 
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[5] Approximate the integral over k modes in equation (14) as a sum to obtain the total anisotropy. Since 
the integrand is oscillatory, unless a large number of k values are taken, spurious oscillations will occur 
in Ci, both in the analytic and numerical calculations. This is in practice not a problem since weighting 
from any realistic experimental window function will automatically smooth out these oscillations. Rather 
than take extra computational time for purely aesthetic reasons, we employ only of order 100 k values 
between the present horizon and the damping scale at recombination but smooth the final results in £ 
for display. One can verify that this is a valid procedure by increasing the number of k divisions. 

[6] Normalize Ci to the COBE DMR detection at large scales. Here we have used the rms 10° value 
{AT/T)Iqo = X;(2^+ l)QW£/47r = 1.25 x 10"!° (Bennett et al. 1994) where the COBE window function 
is Wi = exp[—£(£+l)a'^] with a = 0.0742 as the gaussian width of the 10° beam. For Harrison-Zel'dovich 
CDM models, one may alternatively normalize to Qrms-PS = To(5C2/47r)^''^ = 20yuA' (Gorski, et al. 
1994). 

On the other hand, if only a rough estimate of anisotropics in the CDM model (20% in temperature fluctu- 
ations for scales up to the second Doppler peak) is needed, the following quick and easily coded procedure 
can be used: 

[a] As [2] above, evaluate Oq and ©i in the large scale limit but employ this for the small scales as well. This 
amounts to a 10 — 20% discrepancy in the oscillation behavior in a CDM model where R(ri^,) <^ 1. Use the 
simple damping factor from the instantaneous recombination approximation T>(k) = exp{ — [fc/fc£)(r?H.)]'^} 
with kjy [t]) from equation (B-16). See Fig. 10 for the error this causes. 

[b] Free stream the solution by approximating it with equation (A-23) (10% AT/T errors, see Fig. 7). Follow 
steps [5] and [6] above. 
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Table 1: Commonly Used Symbols 



Symbol 


Definition 


Equation 


At 


Total density fluctuation 


(A-1) 





CMB temperature fluctuation 


(1) 


©0 


CMB Monopole fluctuation 


(1) 


©1 


CMB Dipole fluctuation 


(1) 


e>t 


CMB ftli multipole fluctuation 


(1) 


n 


Anisotropic stress perturbation 


(A-3) 


$ 


Gravitational (Newtonian) potential 


(A-1) 


* 


Gravitational (curvature) potential 


(A-2) 


T] 


Conformal time 


(1) 


m 


Present conformal time 


(11) 


ri* 


Recombination conformal time 


(C-2) 


(XT 


Thomson cross section 


(1) 


r 


Thomson optical depth 


(1) 


V 


Temperature damping factor 


(13) 


A 


Initial power spectrum 


(A-20) 


Ci 


Anisotropy power spectrum 


(14) 


F 


Gravitational driving force 


(5) 


Nt 


Neutrino lib. multipole 


(A-3) 


R 


Normalized scale factor 3pi,/4pj 


(3) 


Ud 


n = decaying mode 


(A-7) 


Ug 


n = growing mode 


(A-6) 


Vt 


Total velocity amplitude 


(A-4) 


Vb 


Baryon velocity amplitude 


(1) 


a 


Scale factor 


(1) 


(IQ 


Present scale factor 


(1) 


Oeq 


Equality scale factor 


(D-8) 


Cs 


Photon-baryon sound speed 


(6) 


u 


Neutrino fraction Pv/ifiv -\- P~i) 


(A-13) 


k 


Fourier mode wavenumber 


(1) 


ko 


Diffusion damping wavenumber 


(10) 


keq 


Equality horizon wavenumber 


(D-8) 


ks 


Solution switching wavenumber 


(B-10) 


i 


Multipole number 


(2) 


He 


Electron number density 


(1) 


r-s 


Sound horizon 


(7) 


Xe 


Electron ionization fraction 


(1) 


overbar 


Undamped (pressureless) solution 


(A-6) 


overdot 


Conformal time derivative 


(1) 


overhat 


Undamped (diffusionless) solution 


(8) 
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